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According to Pruessner and Peters [Phys. Rev. E 73, 025106(R) (2006)], the finite size scaling 
exponents of tlie order parameter in sandpiie modeis depend on the tuning of driving and dissipation 
rates with system size. We point out that the same is not true for avalanches in the slow driving 
limit. 
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In a recent paper Pruessner and Peters investigated 
the relation between self-organized criticality (SOC) in 
sandpiie models and absorbing state phase transitions, 
on the basis of an analogy with standard equilibrium crit- 
ical phenomena, in particular the two-dimensional Ising 
model [1]. According to Ref. [H only a careful choice 
of the system size dependence of the driving and dissi- 
pation parameters would yield the scaling results of the 
underlying phase transition. Here we point out that this 
reasoning does not apply in particular when one studies 
"SOC variables" such as various measures of avalanches 
in the slow driving limit, as is traditionally done in this 
context. This is confirmed by numerical simulations pre- 
sented below. 

In SOC sandpiie models, the steady-state is main- 
tained by a balance of dissipation e and driving h. In 
particular, the control parameter of the absorbing phase 
transition, i.e. the average height of the pile C, evolves 
on the average as 

t^h-ep, (1) 

where p is the density of active sites, the order pa- 
rameter of the absorbing phase transition. In the limit 
— > 0+,e — * 0+, the control parameter flows to the crit- 
ical value Cc and the model shows scale invariance 
The question raised by the authors of Ref. [l[ is how to 
apply finite-size scaling (FSS) to the order parameter p 
if e{L) ^ L~'^ and h{L) ^ are taken to be functions 
of the system size L. To investigate this issue, Pruessner 
and Peters employ a similar "driving" to the Ising model, 
using a fluctuating inverse temperature 6 evolving as 

e = h-(\ml (2) 

where |m| is the absolute value of the magnetization. Us- 
ing this driving, the system size dependence of the order 
parameter, (|m|) ^ L'^^'^, coincides with that expected 
from standard FSS only when lo — k — P/v, where f3 and 
V are the equilibrium exponents for the order parameter 
and the correlation length, respectively. We notice that 



the case cu — k > fi/v corresponds precisely to the slow 
driving limit of SOC, where there is complete time-scale 
separation between driving and dissipation (i.e. for most 
purposes one may even take w — > cx) and wait an inflnite 
amount of time after each driving event). In this limit, 
the effective temperature deflned in Ref. [ij also would 
diverge. The authors, in analogy with the Ising model, 
conclude that the slow driving SOC state should not cor- 
respond to the critical point of an absorbing phase tran- 
sition, in contradiction with the evidence from numerical 
simulations of sandpiie models 

The apparent contradiction disappears when we no- 
tice that avalanche statistics, as typically studied in the 
case of SOC, and the average order parameter studied 
in Ref. Q are not equivalent measures of the criticality 
of the underlying absorbing state phase transition. As 
pointed out by Pruessner and Peters, it is indeed pos- 
sible to tune the system size dependence of the order 
parameter, or average activity, by choosing the size scal- 
ing of driving and dissipation rates appropriately. In the 
slow driving limit relevant for SOC, however, this is just 
a trivial consequence of the drive rate dependence of the 
quiescent periods between avalanches. The avalanches 
themselves are not affected by the drive rate in any way 
as long as it is slow enough such that no new grains are 
added while the system is active. 

In SOC sandpiles, one usually implements open bound- 
ary conditions and infinitely slow drive, corresponding to 
K = 2 and lj ^ oo. The mappings to absorbing phase 
transitions and to dcpinning transitions allows to obtain 
the scaling behavior for any value of k, provided that we 
remain in the time scale separation regime [2[. In gen- 
eral, sandpiie models exhibit FSS forms for the avalanche 
sizes s, of the type 

P{s,i)^s-^^P{s/e)- (3) 

Ts and D are critical exponents related to the underly- 
ing depinning transition and ^(w,/*) is the cut-off scale 
that is determined by the condition of balance between 
dissipation and drive (which also results in the steady 
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state condition p = h/e). One transparent argument to 
compute ^ is to look at the dynamics following the addi- 
tion of a single grain, which gives rise to an avalanche of 
average size (s) dissipating on the average one grain [^. 
Thus one has the condition e(L)(s)(^) — 1 and obtains 

Notice in particular that ^ is not dependent on the drive- 
rate or Lo. 

To confirm this, we consider the Manna sandpile model 
with periodic boundary conditions, slow driving (i.e. 
Lu — oo) and bulk dissipation. The model is studied 
here on a 2d lattice, where for each lattice site i one 
assigns an integer variable Zi (the number of "grains"). 
li Zi > Zc = 1, a "toppling" occurs and the grains are 
redistributed according to Zi —>■ Zi — 2 and z„n = Zn„ + 1, 
where Znn sue two randomly chosen nearest neighbors 
of site i. The dissipation is implemented by removing a 
toppling grain from the system with probability e oc L^". 
Here we consider the two cases with k = 1 and k = 3, 
respectively. In Figure 1 we show that the scaling follows 
the predictions of absorbing phase transitions: the cutoff 
scale of the avalanche size distribution scales according 



to Eq. (lU, regardless of the L-dependence of e. 

To summarize, we would like to point out that the con- 
clusions of Pruessner and Peters are misleading in the 
sense that avalanche statistics in sandpile models indeed 
follows from the underlying absorbing state transition 
whenever one studies the slow driving limit lo — k > (3 /v. 
This is the condition of complete time scale separation 
between driving and avalanche propagation that has been 
recognized already in the early literature as the crucial 
ingredient for SOC in sandpile models 
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FIG. 1: Scaling plots of the avalanche size distributions from 
the two-dimensional Manna model with bulk dissipation. The 
inset presents the collapse in the case n = 1, while in the main 
figure the case k = 3 is shown. Both collapses correspond to 
the value Ts = 1.28. 



